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We studied effects of random potentials and roles of electron-electron interactions in the gapless 
phase of coupled Hubbard chains, using a renormalization group technique. For non-interacting 
electrons, we obtained the localization length proportional to the number of chains, as already 
00 ■ shown in the other approaches. For interacting electrons, the localization length is longer for stronger 

interactions, that is, the interactions counteract the random potentials. Accordingly, the localization 
length is not a simple linear function of the number of chains. This interaction effect is strongest 
when there is only a single chain. We also calculate the effects of interactions and random potentials 
on charge stiffness. 
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I. INTRODUCTION 



One of the attractive topics in mesoscopic one dimensional (ID) electron systems is the disorder effect. While no 
extended state can survive in the presence of random potentials in ID infinite systems, the systems of finite size can 
!h ■ be metallic if the size is smaller than the localization length. It is always important to take a serious look at interplay 
between random potentials and electron-electron interactions when we discuss the transport properties of mesoscopic 
metallic wires, which are now available as an achievement of the recent technological progress. An interesting example 
is a mesoscopic metallic ring, which have attracted much attention since a large persistent current was observed even 
with a modest amount of impurities in the ring jjj. A simple study, only taking an account of the impurity effect, 
failed to explain such large current ||, and it is necessary to consider electron-electron interactions. 

It is in general quite difficult to take correctly account of interactions. One of the exceptions is ID systems, where 
interactions can be treated rather rigorously with help of the bosonization techniques. A renormalization group (RG) 
calculation for ID systems with impurities was performed by Giamarchi and Schulz ||]. One of their interesting 
l— ~ 1 , results is that partilce-particle interactions of spinless Fermions would enhance the disorder effect and help the system 
localize, while the interactions of spinning Fermions {e.g. Hubbard type interactions) would counteract the disorder 
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The RG calculation showed that the effective backward impurity scattering IF of a spinless Fermion system is 
renormalized as dW/dl = (3 — 2K p )W where K p is the Luttinger parameter of the charge mode ||. The larger 
repulsive interactions give rise to the smaller K p and therefore to the stronger effective impurity potentials. This is 
due to the enhancement of charge density wave (CDW) correlations in the ground state of the system. The repulsive 
interactions would enhance the CDW correlations, which make the system easily pinned by the impurities. For 
electrons with spins, however, the renormalized impurity potential is given by dW/dl = (3 — K p — K a — y)W, where 
K u (y — p, a) is the Luttinger parameters of charge and spin modes respectively and y measures the backward 
scattering strength between electrons of opposite spins [00). For Hubbard model with small U, this RG equation 
becomes 
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where vp is the Fermi velocity. Here the electron-electron interactions would screen the impurity potentials. Since a 
spin density wave (SDW) correlation is dominant in the ground state of the repulsive Hubbard model, the interaction 
U makes the particle density uniform, and therefore make the coupling of the density to the impurities weak. 

These opposite roles of interactions in the spinless and spining Fermions was later checked numerically in connection 
with the problem of persistent current (^-0] . Although the suppression of the effective disorder strength due to electron- 
electron interactions in the models with the spin degrees of freedom might explain the observed large persistent current 
in the disordered metal rings, we should carefully consider the fact that the metal rings have the finite cross section 
and therefore have the finite number of channels. 

Ladder systems have recently attracted a wide range of attentions. Disorder effect on electronic ladder systems is 
however not widely investigated so far, and there are few studies which calculated RG equations for two-leg electronic 
ladder with impurity potentials |^-[l^]. If we just focus on the role of interactions in the presence of impurities, it is 
shown that the role changes drastically depending on the phase; the impurity effect is enhanced by electron-electron 
interactions in C1S2 phase (CnSm means that n charge modes and m spin modes are gapless) || while it is suppressed 
in C2S2 phase (with all modes being gapless) [jlO| . Remember that the ground state of a single chain is gapless and the 
impurity effect is suppressed by the interactions. So, as far as we concern the gapless phase, the role of interactions 
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acting against random potentials does not alter even when we increase the number of channels from one to two. But 
the point we should note is that the suppression of the effective disorder strength due to interactions is smaller in two 
channel ladder than single channel chain |10|. If this tendency continues in the systems with more channels, we could 
not expect the large interaction-driven suppression of electron localizations in a thin but finite cross-section wire. It is 
therefore of importance to see in the increased number of channels whether electron-electron interactions counteract 
random potentials and if so how large the effect is. 

In this paper we first show in Sec. 2 how the disorder effect on a non-interacting electronic ladder changes as 
the number of chains is increased. In Sec. 3 we turn on electron-electron interactions and investigate the interplay 
between the interactions and the impurity potentials in one-, two-, and three-chain ladders. Section 4 is devoted for 
summary of the paper. 



II. NONINTERACTING DISORDERED CHAINS 



We start from a tight-binding Hamiltonian of coupled chains with open boundary condition, 

Hq = -t'^ j c { miji c mal+ i -t±^] cl ncrl c m+lai + h.c, 



(2) 



where m and i are the chain and site indices, respectively. This hopping terms are diagonalized by a unitary trans- 
formation c m = '^2 a V ma a a and they form bands. The matrix V for one, two, three-chain systems, for example, is 
given by 
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1 2 V2 
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(N ch = f), 

(N ch = 2), 
(N ch = 3), 



(3) 



where N c h is the number of the chains. Following the standard bosonization procedure, we get the Tomonaga-Luttinger 
type Hamiltonian 



J Z7T i^av 
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where v = p,a represents charge and spin mode respectively, and a is the band index. With no interaction, the 
parameters K are equal to 1 and u are the Fermi velocities of each band. 
The impurity scattering potential term is originally written in the form 
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where V* is the complex conjugate of V . After the bosonization we get the forward and backward scattering terms. 
Since the forward scattering does not affect Hq |9j , we only consider the effect of the backward impurity scattering. 
The backward scattering terms are given by 



H b = — [dxY] ^e-^+^ +s ^ + ^^cos[{(0 Qp - d Pp ) + S (6 aa - 9 Pa )}/V2] + h.c, 

ap,s 



where a is the lattice constant. £ is assumed to be a linear combination of Gaussian random variables and then the 
following general form of the action will appear in the replica trick method by integrating out the random variables, 
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2 



where i, j are the replica indices. In order to discuss the disorder effect of various coupled chains on the same basis, 
is a Gaussian random variable and satisfies < n™?n^ v >= W5 m idij. Hence we have 



z "s =wj2 v: a Vn0v ni v: 5 . (9) 



The RG equation for Z"j? is 
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where we used the fact that JT's are all equal to 1 for free particles. The solution of the RG equation is Z"g(l) = 

e l WJ2 n V n a V npV ni V* s . Note that We 1 = J2p z 2p( 1 )- The RG equation stops when Z"£(l) reaches ~ W 2 /a, and it 
is when I ~ log(L; oc /a) where L; oc is the localization length. Putting all together, we get 

Lloc ~ N ch ^r. (11) 

On the other hand, the electron scattering rate is given by r _1 ~ Wp(ep), where p(ep) (~ Vp 1 ) is the density of 
states at the Fermi level, and hence the mean free path l e = vpr ~ Wf 2 /W . Then we get the known relation ElK2[, 

Lloc ~ N c hl e - 

In order to see the disorder effect on the transport properties, we next calculate charge stiffness D, which measures 
the strength of Drude peak, a — DS(lu) + a reg . Note that D is also a measure of the persistent currents for small flux. 

The external flux couples with the current ^2 as jas, where a is the band index. Since ^2 as jas oc X^Hapi charge 
stiffness is proportional to ^ a K ap u ap . Ignoring the irrelevant numerical factor, we define charge stiffness per channel 
by D = (1/N c h) X] a K ap u a p. For non-interacting and non-impurity ladder, D is the averaged Fermi velocity vp. 

In the presence of impurities the RG equation for the charge stiffness is given by 

dl - l^ u l^N ch \ Z ^ (12) 

IT F 

D has no explicit N c h dependence and the effect of increasing channel number could appear only through Vp 1 , which 
is usually weak. For a finite system of size L, D is given by 

D ~D - -~Vp I WL, (13) 

where Do is a constant. For a given size L, the charge stiffness D is smaller for stronger disorder, and for a given 
disorder W, D is smaller for a larger system. The former is trivial, and the latter is because Li oc /L becomes smaller 
as L increases with W fixed and hence the system will be in more localized regime. This will be clear if we rewrite 
Eq.(||) as 

D ~ D — ~N ch VF~y—' (14) 

We cannot simply extend these discussion to N c h — > oo, that is, to two dimensions (2D), because we constructed 
the RG theory based on the Tomonaga-Luttinger type bosonized Hamiltonian with well-defined subbands, which is 
not an appropriate basis for 2D system. 



III. INTERACTING DISORDERED CHAINS 



In the previous section we showed that the localization length is proportional to the number of chains, and the 
effect of random potentials becomes weaker as the chain number increases. What happens if we turn on interactions 
between the particles. The original Hamiltonian is then Eq.(0) plus the Hubbard-type interaction term, 
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A part of the interaction terms can be combined with the kinetic term to give 

/dx 
— [u ru K ru (irXl rv ) 



2 ' ^{dAr.fl 



(15) 
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where v = p,a represents charge and spin mode respectively, and r is the new band index assigned by the unitary 
transformation (f> aL , — > *^2 r T ar (f) rl , and Ii a u — > ^ r T ar Ii ru , where 
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(N ch = 1), 



(iVch = 2), 



(N ch = 3). 



(17) 



This unitary transformation is necessary to diagonalize Hq. Since a part of interaction terms are included in Hq, the 
parameters K and u have different values from the ones of the noninteracting case. Hereafter we use Greek letters 
a, (3, ■ ■ ■ for the old band index and Roman letters r, ■ ■ ■ for the new band index. 

The other remaining interaction terms are represented by Hi . Since there are still a lot of terms in Hi , we do not 
write them down all here and we only note that the interaction matrix element g a f3-yS in the commonly used notation, 
representing the scattering from (a,/3) bands to (7,<5) bands, is given by 



9af3jS 



IS 



(18) 



in the case of Hubbard model, although some of them, such as g a pf3a and g a pf3p, should be set to zero due to lack 
of the momentum conservation. As stressed in Ref. |l3| ] , one should not forget the sign factors of g which come from 
Majorana Fermion operators. In the present case, the backward scattering interaction g^ga, for example, has to be 
multiplied by —1. 

The presence of the interaction terms made the phase diagram quite rich p3[-pl|, and the interplay between 
interactions and random potentials is strongly phase dependent as stated in Sec. 1. In order to construct a systematic 
view with changing the number of chains, we should pick up a common phase, gapless phase, where no gap opens 
both in the charge and spin modes. There is also another reason to study gapless phase. When the forward impurity 
scatterings are considered, Fujimoto and Kawakami showed |9) that the electron-electron interactions are effectively 
suppressed and as a result the charge and spin gaps collapse. In this situation the remaining task is to study the 
effect of the backward impurity scatterings on the gap-collapsed (gapless) phase. 

In the gapless phase the action presenting the backward impurity scatterings is again given by Eq.(0). The RG 



equation for Z"{? is 



,// E 3 - \ E« ( T - - T ^f + (T ir - T Sr f}{K^ + K-J) 

r 

+{(T ar + T Pr ) 2 + (T rr + T Sr f}(K rp + K ra )}]Z^ 

1_ V n {1) 7 aP 



(19) 



a 1 S3 1 



where g^ represents the backward scatterin g in teractions. Note that, for noninteracting electrons, we have K = 1 
and g = and then Eq. @ reduces to Eq. @. Also, when N ch = 1, Eq. (|l|) reduces to dW/dl = (1 - U/u a ir)W 
which is just what we mentioned in Eq. (|l|). 

Since the renormalization of Z"g has a strong K dependence, the localization length and other physical quantities 
would changes in accordance with K in a complicated way. To make the story simple, we consider the weak interaction 
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limit and use K = 1 as in the noninteracting system. Also we assume the random impurity potentials are so weak 
that we can forget about the renormalization of K . The RG equation ([l9]) can then be written in the following matrix 
form, 

f = (l-0*. (20) 

The element of the N c h x N c h matrix Z is Z^g where i = (a, (3) and j — (j,6). The element of the 

matrix G is given by ^j^5™/3<5 wriere 9^ is now assumed to include the sign created by the Majorana Fermion 
operators mentioned above. The signs of the elements of G are not always plus and it is not in general trivial whether 



the presence of G in Eq. (20) would weaken the growth of Z, namely, whether the interactions would suppress the 
impurity effects. So let us see it in more detail. 

Since we assume the interactions are weak and ignore the renormalization of G, the solution of Eq. (EOT) is 
Z{1) = e(^ G )'Z(0) ~ e\\ - Gl)Z(0). Therefore Tr(Z(l)) ~ e l {N ch W - lTr(GZ(0))} ~ N ch We^- cu ^^ 1 , 
where C = u a irTr(GZ(0))/(N ch WU). Since, when / - log(L loc /a), Tr(Z(l)) ~ N 2 h W 2 /a, we have N ch W 2 /Wa ~ 
(L[ oc /a) 1 ~ cu / VF7T . Then the localization length of the interacting system has the form, 

1+CU/vFn 

(21) 





where L^ c = N c hVF 2 /W is the localization length of the noninteracting system. Because of the presence of U, the 
localization length is no longer simply proportional to N c h- If we forget about the sign of g a p~/8 created by the 
Majorana Fermion algebra, discussed above, C = 1 for any N ch . Taking the sign into account, however, C becomes 

(N ch = 1), 

C = { 0.25 (N ch = 2), (22) 

(Nch = 3). 

Therefore the interaction U always makes the localization length longer, namely, the interactions have a derealization 
effect in the gapless phase of the coupled chain systems. C is not a monotonic function of N c h as far as 1 < N c h < 3. 
and hence it is not always true that the interaction effect becomes weaker as N c h increases. Anyway, since it is easy 
to prove C < 1, the derealization effect is strongest when N c h = 1. 

Next we study the charge stiffness D. As noted in the previous section, the external flux couples to ^2 as j as oc 
J2 a = J2ar ^ar^rp- Recall that the Greek letter a represents the old band index assigned before the operation of 
T, and the Roman letter r is the new band index. Since T ar — 1 (when N c h = 1), \/25 r i (when N c h = 2), y/3S r 2 
(when N c h — 3), the charge stiffness is given by D = u rp K rp where r — 1 (when N c h = 1, 2), r = 2 (when N c h = 3) 
with the irrelevant numerical factor being omitted. The RG equation is then given by 

dD _ 2au 2 rp ^ - 

dT--^T2 ? ^ r - (23) 

7T 

where r = 1 (when N c h — 1,2), r — 2 (when N c h = 3). In the second line of Eq. ( p3| ) the interactions are assumed to 
be small. For a finite system of size L, 

D ^ Do - f^ W ~i W (L\ . (24) 

Again, the interaction U suppresses the random potential effect and gives larger value of the charge stiffness that in 
the noninteracting case. 



IV. SUMMARY 



We investigated the role of electron-electron interactions in the coupled Hubbard chains with random potentials. For 
noninteracting systems, a RG calculation shows that the effective strength of the impurity potentials grows towards 
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the strong coupling limit and the localization length is proportional to the number of chains N c h- In the presence of 
interactions between the particles, the role of the interactions changes from the phase to phase as previously shown 
in Ref. [p|JlO|], and therefore we only focused on the gapless phase of N c h chain systems. Based on the RG calculation, 
we showed that the interactions always reduce the effective strength of the impurity potentials. The degree of the 
reduction is sensitive to N c h, and the counter-effect of the interactions against the random potentials is strongest 
when N c h = !• 
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